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We search for dynamical magnetoelectric phenomena in three-dimensional correlated systems with spin- 
orbit coupling. We focus on the antiferromagnetic insulator phases where the dynamical axion field is realized 
by the fluctuation of the antiferromagnetic order parameter. It is shown that the dynamical chiral magnetic 
effect, an alternating current generation by magnetic fields, emerges due to the time dependences of the order 
parameter such as antiferromagnetic resonance. It is also shown that the anomalous Hall effect arises due to the 
spatial variations of the order parameter such as antiferromagnetic domain walls. Our study indicates that spin 
excitations in antiferromagnetic insulators with spin-orbit coupling can result in nontrivial charge responses. 
Moreover, observing the chiral magnetic effect and anomalous Hall effect in our system is equivalent to detecting 
the dynamical axion field in condensed matter. 
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Introduction .— Antiferromagnets (AFMs) have attracted 
much attention from the viewpoints of both purely scientific 
and applied research. Theoretically, it is known that antifer¬ 
romagnetic (AF) phases are often favored in systems with 
strong on-site repulsive interactions. In the vicinity of AF 
phases, the emergence of exotic phases and phenomena such 
as high-temperature superconductivity and spin liquids is well 
acknowledged [1-3]. On the other hand, AFMs have recently 
been studied intensively in the field of spintronics [4-7], as a 
possible new class of materials alternative to ferromagnets. 
These studies suggest that the staggered magnetization can 
play essential roles, although AFMs had not been considered 
suitable for practical use due to the lack of net magnetization 
unlike ferromagnets. 

Recent intensive and extensive studies have revealed the 
importance of spin-orbit coupling (SOC) in condensed mat¬ 
ter. Especially, the discovery that strong SOC is essential to 
realize topologically nontrivial phases opened a new direction 
in modern physics [8-10]. Since topological invariants are de¬ 
termined from electronic band structures, studies of topologi¬ 
cal phases started as a single-particle problem. Subsequently 
many-body effects in topological phases and spin-orbit cou¬ 
pled systems have become an attractive subject [ 11 , 12 ], and 
the emergence of novel phases such as the topological Mott 
insulator [13] (or more generally the fractionalized topolog¬ 
ical insulators [14]) and the Weyl semimetal [15] has been 
predicted. As for novel phenomena in spin-orbit coupled and 
correlated systems, for example, the axionic polariton, a total 
reflection phenomenon of light, has been suggested [16]. 

In this Letter, we study electromagnetic responses of an¬ 
tiferromagnetic insulator (AFI) phases in three-dimensional 
(3D) correlated systems with SOC. We explore dynamical 
magnetoelectric phenomena where the staggered magnetiza¬ 
tion plays essential roles. We show that, in the presence of 
SOC, spin excitations in AFIs can result in nontrivial charge 
responses, as a consequence of the realization of the dynam¬ 
ical axion field. First we show the emergence of the chiral 
magnetic effect (CME), an electric current generation by mag¬ 


netic fields [17], in the AEI phase. The CME was originally 
proposed in gapless Dirac fermion systems [17], and its pos¬ 
sibility has been discussed in Weyl semimetals [18-23]. In 
contrast to preceding works, we propose the dynamical real¬ 
ization of the CME in gapped systems. We also show the 
occurrence of the anomalous Hall effect (AHE) in the AEI 
phase. It is known that the AHE occurs usually in ferromag¬ 
netic metals [24], while the AHE arising from nontrivial spin 
textures has been studied in frustrated or noncollinear AEMs 
[25, 26]. We propose that spatial variations of the staggered 
magnetization lead to the AHE. 

Realization of the Dynamical Axion Field and its 
Consequences .— Let us consider 3D electron systems having 
both on-site interactions and SOC, such as 5d transition metal 
oxides [12, 13, 27]. We focus on systems that become mag¬ 
netically ordered Mott insulators when on-site interactions are 
strong, while they are topological band insulators when on¬ 
site interactions are weak. Once a magnetic order is formed, 
the mean-field approximation of the interaction term can cap¬ 
ture the essential physics of the system. In this work, we 
particularly consider AEIs whose mean-field lattice Hamil¬ 
tonian is given by ‘7T(fc) = eo(fc)l ■+■ Rfi(k)a^. Here 
k - {k\,k 2 ,kf) is a wave vector in the Brillouin zone, 1 is 
the 4x4 identity matrix, and the 4x4 matrices satisfy the 
Clifford algebra {a^,ay] - 26^v with as = aia 2 asa 4 . The 
Hamiltonian of this form can be realized, for example, in the 
AEI phases of 612803 family doped with magnetic impurities 
such as Ee [16] and transition metal oxides with corundum 
structure such as a-Ee 203 [28]. In this case, we can derive 3D 
massive Dirac Hamiltonians of the form 

’kfefiiq) - qioi H- q2a2 + q^aj, H- -H Ms/as (1) 

around some momentum points Xf, where q - k - Xf. The 
subscript / indicates the valley degrees of freedom. The 
kinetic term q/ja/j is spin-dependent as a consequence 
of SOC. Moa 4 is a mass term with time-reversal and parity 
(spatial inversion) symmetries induced by SOC, and MsfUs 
is a mass term with broken time-reversal and parity symme- 
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tries induced by mean-field AF order parameter. We require 
that the system is a topological insulator when Mq > 0 and 

Mif = 0 . 

In what follows, we consider consequences arising from the 
existence of the Ms/as mass term. The effective action of the 
system in the presence of an external electromagnetic poten¬ 
tial is written as 

5efF = r dtd^r'^ ^f{r, t) t), (2) 

where t is real time, (/^/(r, t) is a four-component spinor, ij/f - 
Df, ^ + ieAf,, M'^ = + (Ms/)^, cos Of = 

MqIM'j., sinfly^ = -MsflM'j-, and we have used the fact that 
a 4 = y“, as - and a,- = y^y-i (J =1,2,3). By applying 

the Fujikawa’s method [29] to the action (2), the theta term is 
obtained as [30] 

5e= r dtcPr^OE- B, (3) 

J znh 

where 9 - |[1 -i- sgn(Mo)] - tan“'(M 5 y/Mo), and E 
(B) is an external electric (magnetic) field. From this ac¬ 
tion, we obtain the magnetoelectric responses expressed by 
P - 9e^I{2nh)B and M - 6e^ l{27Th)E, with P the electric 
polarization and M the magnetization. In 3D time-reversal 
invariant topological (normal) insulators, 9 - n {6 - 0) [36]. 
However, the value of 9 can be arbitrary when time-reversal 
and parity symmetries of the system are broken [37-39]. Fur¬ 
thermore, when the value of 9 depends on space and time, 
it can be said that the dynamical axion field is realized in 
condensed matter [16]. Some consequences of the realization 
have been studied so far [16, 40]. 

Notice that, when the dynamical axion field is realized, the 
theta term can be rewritten in the Chern-Simons form as 

/ 2 

dtd^r^^e^''P\df,9(r,t)]AydpAA. (4) 

Then the induced four-current density f can be obtained from 
the variation of the above action with respect to the four- 
potential ^ = -|^[5;;6>(r,f)]e''''^'*5pA4. The in¬ 

duced current density is given by [41, 42] 

r. 

j{r, t) = — \9{r, t)B + V6»(r, t)xEj, (5) 

where 9 — d9(r, f)/dt. The magnetic-field-induced term is the 
CME [17]. The electric-field-induced term is the AHE, since 
it is perpendicular to the electric field. The induced current 
of the form (5) has been also studied in Weyl semimetals [18- 
21 ], where the chemical potential difference between the band 
touching points and the separation of the points in momentum 
space are required for the CME and AHE, respectively. How¬ 
ever, the existence of the CME in Weyl semimetals is still be¬ 
ing discussed theoretically [18-23]. Note that the situation we 
consider in this paper is completely different, since the system 


is gapped, i.e., the above conditions required in the case of 
Weyl semimetals are not needed. 

Theoretical Model .— To study the induced current (5) more 
concretely, let us consider a 3D lattice model with SOC and 
electron correlations. The model we adopt is the Eu-Kane- 
Mele-Hubbard model on a diamond lattice at half-filling, 
whose Hamiltonian is given by [43-45] 

^ + 'll Z ^ ^ 

(ij).cr ((i.j)} 

i 

where is an electron creation operator at a site i with spin 
i)j «!o- = Cicr, and a is the lattice constant of the fee 
lattice. The first through third terms represent the nearest- 
neighbor hopping, the next-nearest-neighbor SOC, and the 
on-site electron-electron interaction, respectively, and 
are the two vectors which connect two sites i and j on the 
same sublattice. Namely they are given by two of the four 
nearest-neighbor bond vectors, cr = (ctj, 0 - 2 , (T 3 ) are the Pauli 
matrices for the spin degrees of freedom. We introduce a lat¬ 
tice distortion such that = t + 6t\ for the [111] direction and 
tjj - t for the other three directions, which induces a bandgap 
of 2\Mq\ {Mq = dti) in the noninteracting spectrum. 

We perform the mean-field approximation to 
the interaction term as Nu = 17 2, niT«/i ~ 

t7 2/ [<n,i>n,T + («iT>”4 “ 

xcj^c;| -I- SOC breaks spin SU(2) symmetry 

and the orientations of the spins are coupled to the lattice 
structure. Hence we should parametrize the AE ordering 
between the two sublattices A and B in terms of the spherical 
coordinate (n, 9, ip): 

{Sr a) = -{Si'b) - {n sin 9cos tp, n sin 9 sin (p,n cos 9) 

= Tiie^ -H n2ey + {= n), 

where <5'/-^) = \{cl^^(TapCi’f,p) {pi ^ A, B) with i’ denoting 
the i'-th unit cell. In the following we consider the ground 
state given by {hq, 9q, ipo). The low-energy effective Hamil¬ 
tonian of the AEI phase is written in the form (1); //eff = 
IiqIi/=i.2,3'A}q'77/(q)iA/<j- where tj/fg is a four-component 
spinor [45]. Therefore the value of 9 in the AEI phase of the 
Eu-Kane-Mele-Hubbard model is given by [45] 

= ?[1 + sgn(Mo)] - ^ tan^' (t/n//Mo) . 

/=1,2,3 

Erom this equation, we see that the dynamical axion field is 
realized by the fluctuation of the AE order parameter nf, i.e., 
by the spin excitations [16, 28]. 

Dynamical Chiral Magnetic Effect .— Eirst we focus on the 
magnetic-field-induced term in Eq. (5), i.e., the CME in the 
AEI phase; 

, s V UMq . / . s,TJ 
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FIG. 1. (Color online) (a) Schematic figure of the dynamics of mx 
and mg in the AF resonance state, (b) A possible experimental setup 
to observe the CME in our system. A static magnetic field B is 
applied along the easy axis of the AF order, = no/|no|. When the 
external microwave (i.e., ac magnetic field) frequency cj is equal to 
the resonance frequencies tu+, the induced alternating current JcmeCO 
will be observed. 


Let us consider a case where a microwave (i.e., ac magnetic 
field) is irradiated and a static magnetic field B - Z?e„„ is ap¬ 
plied along the easy axis of the AF order. Here e„„ = no/|no| 
is the unit vector parallel to the easy axis. The dynamics of 
the sublattice magnetizations (5 '/'a) = rriA and {Si's) = 'nriB 
can be described by [46, 47] 

niA^mAX {-ujmB -i- [g/tsB -H ajAiniA ■ e„J] e„„), 

rfiB = tub X {-cojmA -b [subB + coAirriB ■ e„J] e„„}, 

where atj and oja are the exchange field and anisotropy field, 
respectively. We write niA and tub as niA - noe„^+6mA±e“^‘ 
and rriB - -noe„o + dniB^e"^'. Up to linear order in 6mA{B)± 
(i.e., \6mA(B)±\ 1), we obtain the resonance frequencies as 

u) - io± - gusB + V(2wy -H UA)tjOA- In the resonance state, 
where all the spins are precessing around the easy axis with 
the same frequency (or the AF order parameter is 
described as 

n±(f) = [niAit) - mB(f)]/2 noe„„ -H (11) 


From the relation such that n - rii e^ -i- n 2 ey + nj,e^, we have 
Ml = nosin0ocos^o + dncosa»fcos0ocos^o~<5nsina»f sin^Oj 
n 2 - no sin 9q sin + 6n cos cot cos Oq sin ipQ + 6n sin cot cos (fQ, 
and n3 = no cos Oq - 6n cos lot sin Oq. Substituting these quan¬ 
tities into Eq. (9), we obtain the analytical expression for 
_7 cme( 0- Especially, in the vicinity of the phase boundary 
where UnjIMo 1 [49], Eq. (9) is simplified as 


Jcme(0 - 


Inh Mo ^ 

“ a=± 


tOadria sin (toj 


a), 


( 12 ) 


where Di = ^Jp^ + and tana = qjp with p - (cos(,oo + 
siny?o)cos0o - sin 00 and q - siny?o - cosi/Jo- Equation (12) 
means that an alternating current is induced by the AE res¬ 
onance. Schematic figure of a possible experimental setup 
to observe the CME in our system is shown in Eig. 1(b). 

is a function of the external microwave frequency co with 
Lorentzian structure, i.e., dn±(co) ~ a/[(co - co±)^ + a^] with 
a being a constant. Therefore two peaks will appear in the 
intensity |jcME(w)l- 

Here let us estimate the maximum value of the CME (12): 
IjcMElmax = Substituting possible values 

Uno/\Mo\ ~ 0.1 [50], |Z)i| ~ 1, (5n±/no ~ 0.02, and 
co± ~ 500 GHz at B ~ 1 T [51, 52], we have IjcMElmax ~ 
1 X 10"^ A/vn?. This value is experimentally observable. It 
should be noted that the current is adiabatically induced in 
the gapped phase as in the case of the quantum Hall effect or 
the topological charge pumping effect [53]. Hence there is no 
energy dissipation, unlike the conventional transport regime 
which causes the Joule heat. 

The CME was originally proposed in the ground states of 
massless Dirac fermion systems as a direct current generation 
by static magnetic fields due to the presence of the chemi¬ 
cal potential difference between band touching points [17]. If 
such a static CME exists in realistic materials, there will be 
substantial potentials for its applications, since the current is 
dissipationless. However, the existence of the CME remains 
a theoretically controversial subject in Weyl semimetals [18- 
23]. As discussed in Ref. 19, the possibility of the static CME 
would be ruled out in crystalline solids (i.e., lattice systems), 
regardless of the presence or absence of energy gaps. In con¬ 
trast to preceding works, our study proposes that the CME oc¬ 
curs dynamically in insulating systems, which requires time 
dependences of the AE order parameter nj caused by external 
forces such as AE resonance state. 

Anomalous Hall Effect .— Next we focus on the electric- 
field-induced term in Eq. (5), i.e., the AHE in the AEI phase: 


Here we have neglected the difference between the angles 
Oa - {\6mAiMno) and Ob — tan“*(|dmBy|/«o)- The ratio 
OaIOb is obtained as OaIOb ~ (1 -h ^IuaIcoj)^ [46]. Typically, 
the ratio <jJaI<jOj{~ KfJ), with K and J being the strength of 
the anisotropy and exchange coupling, respectively, is of the 
order of 10“^ to 10“^ [48]. Therefore we see that a; -tub 
and thus rriA ~ -rriB- Illustration of the dynamics of ttla and 
tub in the AE resonance state is shown in Eig. 1(a). 


UMo 

JAHe(t’, t)^-— > TA-^ ^ (13) 

iTih ^ Mq+ (Unjff ' ^ 

In order to obtain a concrete expression for this curTent, let us 
move on to a new Cartesian coordinate (X, Y,Z). We consider 
a one-dimensional AE texture of length L along the Z direc¬ 
tion, an orientational domain wall [54, 55]. As shown in Eig. 
2(a), the AE order parameter n(r) = [mA{r) - mB('r)]/2 
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FIG. 2. (Color online) (a) Schematic figure of a one-dimensional 
AF texture, an orientational domain wall of length L. The AF order 
parameter n(r) = [mA(r) - rngfr)]/! at the two edges has a relative 
angle 5. (b) A possible experimental setup to observe the AHE in our 
system. A static electric filed E is applied perpendicular to the AF 
order. There is a relative angle 6 in the magnetization directions of 
ferromagnet 1 (FMl) and ferromagnet 2 (FM2). 

at the two edges has a relative angle 6. Namely, we have 
B{Z - 0) - 00 and 0(Z — L) - flo + (^ in the original 
spherical coordinate. A static electric filed E is applied per¬ 
pendicular to the AF order as E - Eyey- For simplic¬ 
ity, we assume that the system lies near the phase boundary 
where UnfjMo ^ 1 [49]. Noting that only the X compo¬ 
nent survives, we see that Eq. (13) is simplified to be 
^'ahe(^) = 2/ dnfiZ)ldZ. The total current in the X 

direction is given by 

JIhe = dZjl^^iZ) = — -j^Ey, (14) 

where D 2 = Ij f dn f = 2;/[«/(6'o + (5) - «/(0o)] = 

nol V2 sin(^o + |)[sin(0o + d) - sin 0o] + cos(6>o -H d) - cos 0o)- 
The Hall conductivity is estimated as ctxy — ~ 

1 X 10“^ e^/h, since Uno/\Mo\ ~ 0.1 [50] and |Z)2l/no ~ 1- 
Schematic figure of a possible experimental setup to observe 
the AHE in our system is shown in Eig. 2(b). Two ferromag- 
nets with a relative angle 6 in the magnetization directions are 
attached to the AEI [55]. In experiments, the 6 dependence of 
the Hall conductivity will be a direct evidence for the obser¬ 
vation of the axion field. Note that, in contrast to preceding 
works on the AHE in AEMs [25, 26], the AHE studied here 
does not occur in uniform ground states. Namely, spatial vari¬ 
ations of the AE order parameter rif need to be realized by 
external forces. 

Discussions and Summary .— Let us discuss briefly the real¬ 
ization of our predictions in realistic correlated systems with 
SOC. It has been suggested that the dynamical axion field can 


be realized by spin excitations in the AEI phases of 612803 
family doped with magnetic impurities [16] and transition 
metal oxides with corundum structure such as Q;-Ee203 [28]. 
In the same manner as above, we can derive similar expres¬ 
sions for the CME and AHE in these systems. What about 
the possibility in other systems? Eirst of all, time-reversal and 
inversion symmetries of the system must be broken to induce 
the deviation of 6 from 0 or n. Theoretically, the value of 0 can 
be calculated numerically in any insulating systems [37-39]. 
The point is that the emergence of the CME and AHE depends 
on whether 0 is a function of physical quantities such as AE 
order parameter, as in our case. If 0 is a function of a physical 
quantity, then the fluctuation of the physical quantity realizes 
the dynamical axion field. It should be noted that, even if the 
value of 0 is zero in ground states, the realization of dynamical 
axion fields is possible. 

In summary, we have studied theoretically 3D AEIs with 
SOC, focusing on a role of the staggered magnetization. We 
have revealed that, in the presence of SOC, spin excitations 
in AEIs can result in nontrivial charge responses. It is shown 
that the dynamical CME, an alternating current generation by 
magnetic fields, emerges due to the time dependences of the 
AE order parameter. It is also shown that the AHE arises 
due to the spatial variations of the order parameter. These 
two phenomena are the consequences of the realization of the 
dynamical axion field in the AEI phase. The magnetic-field- 
induced and electric-field-induced currents in this study are 
understood as a polarization current in the bulk and a mag¬ 
netization current in the bulk, respectively, which can flow 
in insulators. Observing these phenomena is equivalent to de¬ 
tecting the dynamical axion field in condensed matter. In other 
words, we propose a new way to detect the dynamical axion 
field. 
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Supplemental Material 


1. Alternative Derivation of Eq. (3) 

In this section, we derive the effective action consisting of the Neel field n and an external electromagnetic potential A^. For 
this purpose, it is convenient to adopt a perturbative method rather than the Fujikawa’s method. We start with the total action of 
an antiferromagnetic (AF) insulator described by Eq. (1) in the presence of an external electromagnetic potential 


(SI) 
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By integrating out the fermionic field, we obtain the effective action for n and as 


Z[n,A,] = f [^, = exp] ^ Trln [Go)(1 + Gq/V/)] U exp ^ Tr(lnGo;) + ^ Z -Tr(^o/^/} 

^ / L / / «=i " 


(S2) 


= giWcdnAA 


where Go/ = {iy^d^ - Mq) ' is the Green’s function of the noninteracting part, V/ = -ey^A^ + iy^M^f is the perturbation term, 
and we have used iy^D^ - Mq + iy^M^j - G^j + Vf. Next we use the following identities for the traces of gamma matrices; 

tr(y) = tr(7^) = 0, tr(7^y’') = 4g^'', tJ'(y^y''y^) - 0, tr(7^7''y’7°'7^) = (S3) 

The relevant terms up to the leading order read 

iW^s[n,A^] = ]^'^i:x{GQfiy^ Msff + Z Tr [(-Go/e/A^)^ (Go/iy^Ms/) 


/ 


/ 


(S4) 


The first and second terms correspond to a bubble-type diagram and a triangle-type digram, respectively. 

For concreteness, we consider the case of Msj = Urii, Msg - Un 2 , and M53 = Utis (U is the on-site interaction strength), 
which is applied to the Fu-Kane-Mele-Hubbard model. Flere the Neel field is given hy n - riie^ + n 2 ey + n^e^. The first term in 
Eq. (S4) is given explicitly by 


Tr 


/ .5 _ 2 C d'^q C d'^k ir\i{'fk^+MQ)iy^nf{q)i[y''{k + q)y + MQ\iy^nf{-q)\ 

[Gofiy Msf) - U J ^ (^2 _ _ ^2] 

r d*q r 

J (27r)4j ( 


= 4G-' 


= / 


(k^-M^)[(k + q)2 
d'^k [k^{k -t- qY + Ml]nf(q)nf{-q) 
^ (^2 _ Ml)[(k + qY - Ml] 


(S5) 


d*q 


I(q)nf(q)nf(-q), 


where k^ - g'^''kf^ky — k^k^^ — k^ - k^. We have used Gof(k) — i^y^k^j + Mq)I(JY - Mq), {'/',y^) = 0, and {y'^,y''} - 2^^”. 
Furthermore we introduce the Feynman parameter to combine the denominator as 




da- 


(^2 - Ml)[(k + q)^ - Ml] Jo [12 + a(l- a)q2 - mI^ ’ 


(S6) 


where I - k + aq. On the other hand, in terms of /, we obtain k/^ik + qY - fi - a(l - a)q2 H- (1 - 2a)lftq^. Then the integral I(q) 
can be represented as 


I(q)/i = 4G- 


r d% r' 
J (2;r)4Jo 


da 


-I- a(l - a)q2 - Ml 
[ll + Ml - a(l - a)q2]2 


= A + Bq^+ 0(q‘^), 


(S7) 


where we have Wick-rotated as 1^ = if. The term that contains l^q^^ vanishes, because it is an odd function of 1^. Finally, by 
substituting Eq. (S7) into Eq. (S5), we arrive at the action of the form 


ZTr(Go/!yM5/)^ = - J" dtd^r [(5^n) ■ (B^n) + m^n^j, 


(S8) 


where Ijg - B and rrY - AjB. This action is nothing but the action of the Neel field (i.e., the nonlinear sigma model) [31]. In 
the present low-energy effective model [Eq. (SI)], the information on the anisotropy of the Neel field is not included. On the 
other hand, many (actual) AE insulators have the easy-axis anisotropy. Hence the term nYrY will be replaced by a term like 
m^(n ■ ejyY with denoting the easy axis. 

The second term in Eq. (S4) is the so-called triangle anomaly, which gives the theta term. The final result is [32, 33] 


ZTr [-GofcY^Aijf {Cofir’Msf) =/ f dtd^r^ ^[1 + sgn(Mo)] - Z 

-i I dtd^r— — 6 {r,t)E ■ B, 

J 2nh 


Unf(r,t) 


Mq 


Y''P^d^AydpA2 


(S9) 
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where 0(r, t) — |[1 + sgn(Mo)] - 2/ Unf{r, t)IMo. Note that this value of 0 is consistent with that obtained by the Fujiawa’s 
method 0^ - |[1 + sgn(Mo)] - 2 /t™ since tan“'(x) = x- x^/3 + ■ ■ ■. 


2. From the Fu-Kane-Mele-Hubbard Model to the Heisenberg Model 


In this section, we consider the validity of the dynamics of the sublattice magnetizations described by Eq. (10) in the Fu-Kane- 
Mele-Hubbard model. A mean-field study has shown that the AF insulator phase of the Fu-Kane-Mele-Hubbard model develops 
when t//f ~ 4 (with A/t ~ 0.3) [34]. Hence we may consider a strong coupling picture. When U '3> t and t/ » 4, we can derive 
the Heisenberg model by treating the nearest-neighbor (NN) hopping term H, = Tjii.fi.a Ujc]^Cjcr and the next-nearest-neighbor 
(NNN) spin-orbit coupling term - ‘AiAla^Yi((i,j))‘^\(^ ' (<^,7 ^ djj)Cj as perturbation terms. After a standard second-order 
perturbation calculation, we obtain the NN exchange interaction 


H 


NN 

exch. 


EE 

(ij) o'.o-' 






0-U 


4f?. 

(i,i) 


+ const. 


(SIO) 


Similarly, we can obtain the NNN exchange interaction term. For example, in the case of d}.. - |(1,1,1) and d^. — f (1,1,-1), 
we have cjcr ■ (dj. x df.)cj - |-cj[-cr^: + cry]cj. Then a second-order permrbation reads 


H 


{ 5 iAcj [- 0 -X -H o-y\Cj^ (I iAc\ [-cr^, -I- o-y]ci] 24 


exch. 


o-u 




(Sll) 


We can do the same procedure as above for the other NNN directions. 
Combining these things, the exchange interaction term is written as 


(i,j) ((',]}) a,b 


(S12) 


where = 4-tfjlU and A"? = O(A^IU) (a,b — x,y,z). The fact that A"? contains positive values indicates that the NNN 
exchange interaction can favor AF alignments between the same sublattices. Namely, the NNN exchange interaction competes 
with the NN exchange interaction which favors ferromagnetic alignments between the same sublattices. Such a competition is 
consistent with a numerical result from the weak coupling that shows the increase of the critical strength Uc for the AF ordering 
in the presence of spin-orbit coupling [34]. Furthermore, spin-orbit coupling breaks spin SU(2) symmetry, i.e., the orientations 
of spins are coupled to the lattice structure. More generally, as well known, the presence of spin-orbit coupling leads to magnetic 
anisotropy [35]. Therefore, the effects of the second term in Eq. (S12) can be understood as to modify the value of the NN 
coupling Jij and to determine the ground-state direction of spins. 

In many antiferromagnets, there exists the easy-axis anisotropy such that i/aniso. = -K 2,(5', ■ e„„)^ with denoting the 
ground state direction. For example, in the case of the AF insulator phase of Bi 2 Se 3 doped with magnetic impurities such as Fe, 
the easy axis is the z direction (perpendicular to the quintuple layers) [16]. Then, in the presence of an external magnetic field 
B, the spin Hamiltonian of the system could be written effectively as 

^^eff = Si- Sj-K ^(Si ■ e,p Si, 

(',» ‘ i 


where J includes the effect of spin-orbit coupling and we have assumed isotropic exchange interaction for simplicity. It should 
be noted that the continuum action of this Hamiltonian is given by Eq. (S 8 ) (this is the famous Haldane’s mapping) [31]. 

The spin dynamics is determined from the equation of motion: 

Si^i[H,f,,Si]^SiXfi, (S14) 


where fi is the effective magnetic field given by 

fi = -dH^sIdSi = Sj + 2K(Si ■ e„„)e„„ -H gfisB. 

(iJ) 

In order to consider the AF resonance state, where all the spins are precessing around the easy axis with the same frequency, we 
can replace the spins S'; by the mean-field values itia and mg with A, B denoting two sublattices of a diamond lattice. Finally, 




we arrive at Eq. (10): 
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rhA^niAX {-cjjniB + [gf^sB + coAirriA ■ e„„)] e„„), 

ThB^niBX {-cjjniA + [gf^BB + coAirriB ■ e„„)] e„J , 

where atj - /znn (^nn is the number of NN bonds), cja - 2K, and a static magnetic field is applied along as B - Be„^. Note 
that the easy-axis anisotropy term in Eq. (S13) is not essential to cause the dynamics. What is essential is an external magnetic 
field, as is understood from the resonance frequency to - <jj± - g^iBB + yl{2ioj + a)A)coA- 




